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Abstract 

Starting from the Heisenberg-Euler effective Lagrangian, we determine the photon current and 
photon polarization tensor in inhomogeneous, slowly varying electromagnetic fields. To this end, we 
consider background field configurations varying in both space and time, paying special attention 
to the tensor structure. As a main result, we obtain compact analytical expressions for the photon 
polarization tensor in realistic Gaussian laser pulses, as generated in the focal spots of high-intensity 
lasers. These expressions are of utmost importance for the investigation of quantum vacuum 
nonlinearities in realistic high-intensity laser experiments. 
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I. INTRODUCTION 


As has been recognized already in the early days of quantum electrodynamics (QED), 
the fluctuations of virtual charged particles in the vacuum give rise to nonlinear, effective 
couplings between electromagnetic fields However, so far the pure electromagnetic 

nonlinearity of the quantum vacuum though subject to high-energy experiments { 4 ] has not 
been directly verified on macroscopic scales. The advent and planning of high-intensity laser 
facilities of the petawatt class has triggered a huge interest in proposals to probe quantum 
vacuum nonlinearities in realistic all-optical experimental set-ups; for recent reviews, see 


B 


9|. Prominent examples are proposals intended to verify vacuum birefringence 10ljl3| 


with the aid of high-intensity laser s 114] (c f. also [15|), so far searched for in experiments 


using macroscopic magnetic fields 


Complementary suggestions promote the use 


of time-varying fields and high-precision interferometry 


18 


20]. Other commonly studied 

nonlinear vacuum effects encompass direct light-by-light scattering jl, 21], photon splitting 


island spontaneous vacuum decay in terms of Schwinger pair-production in electric fields 

on optical signatures of quantum 


22 


23]. Besides, recent theoretical studies have focusec 


vacuum nonlinearities based on interference effects 24j 26| . photon-photon scattering in the 
form of laser-pulse collisions |27|, quantum reflection |28|, photon merging [29j, and harmonic 


generation from laser-driven vacuum |30|, |31] 


A central object in the study of such effects is the photon polarization tensor, which 
is known analytically in several limits, namely for homogeneous electric and/or magnetic 



also la 


41| and references therein), and generic plane wave backgrounds 


44] for a more recent derivation and an alternative representation, and 29] 


for a novel systematic expansion especially suited for all-optical experimental scenarios. 
Numerical results for inhomogeneous magnetic backgrounds are available from worldline 
Monte Carlo simulations 


45|, and first analytical insights into the photon polarization tensor 


in various inhomogeneous background held configurations that may locally be approximated 
by a constant were obtained in 28]]. However, the latter approach is manifestly limited to 
certain photon polarization modes as the incorporation of the inhomogeneous held prohle 
into the constant held polarization tensor generically induced violations of the Ward identity. 

In this paper we study probe photon propagation in generic inhomogeneous electromag¬ 
netic helds varying in both space and time. In contrast to previous studies which implicitly 
assumed constant background helds when extracting the photon polarization tensor from 
the Heisenberg-Euler Lagrangian (cf., e.g., [ 12 ], 14]), we treat the background held as man¬ 
ifestly inhomogeneous from the outset. Our main interest is on analytical insights beyond 
the explicitly known held configurations, i.e., beyond constant and homogeneous and plane 
wave background helds. 


2 




















































II. LOW-ENERGY EFFECTIVE THEORY OF PHOTON PROPAGATION 


At one-loop order, the effective Lagrangian in constant external electromagnetic fields 
(“Heisenberg-Euler effective Lagrangian”) |2|, describing the effective nonlinear interactions 
between electromagnetic fields mediated by electron-positron fluctuations in the vacuum, 
can be represented concisely in terms of the following propertime integral 23] (cf. also 


0 


4 ». 


C = 


a 

2tt 


ds 


abcoth(as) cot (bs) — 


-b 2 


( 1 ) 


with the prescription m 2 —>■ m 2 — z0 + , and the propertime integration contour assumed to 
lie slightly below the real positive s axis. Here, m is the electron mass, e is the elementary 
charge, a = e 2 /(47r) is the fine-structure constant, and a = (VJ 72 + G 2 — F) 1 ^ 2 and b = 
(■v/J-" 2 + Q 2 + J 7 ) 1 / 2 are the secular invariants made up of the gauge and Lorentz invariants 
of the electromagnetic field: T = \F IW F' JL ' = \(B 2 — E 2 ) and Q = \F IIV *F IIU = —E ■ B, with 
*F flu _ denoting the dual held strength tensor; is the totally antisymmetric 

tensor, fulfilling e 0123 = 1. Our metric convention is g /u , = diag(—1,+1,+1,+1), and we 
use units where c = h = 1. To keep notations compact we moreover employ the shorthand 
notations f x = f d 4 x and f k = f for the integrations over position and momentum 
space, respectively. 

The effective Lagrangian © is a gauge and Lorentz invariant quantity. Clearly, for 
inhomogeneous background fields additional gauge and Lorentz invariant building blocks 
become available. For slowly varying fields the deviations from the constant held limit can 
be accounted for with derivative terms ~ d a F^. If the typical frequency/momentum scale 
of variation of the inhomogeneous background held is v, derivatives effectively translate into 
multiplications with v to be rendered dimensionless by the electron mass m. Thus, Eq. © is 
also applicable for slowly varying inhomogeneous helds fulfilling — -Cl, or in other words 
- for inhomogeneities whose typical spatial (temporal) scales of variation are much larger 
than the Compton wavelength (time) ~ A of the virtual charged particles. The electron 
Compton wavelength is A c = 3.86 • 10 _13 m and the Compton time is r c = 1.29 • 10~ 21 s. In 
turn, many electromagnetic helds available in the laboratory, e.g., the electromagnetic held 
pulses generated by optical high-intensity lasers, featuring wavelengths of 0(n m) and pulse 
durations of O(fs), are compatible with this requirement. The effective Lagrangian is a 
scalar quantity, and scalar quantities made up of combinations of F^ v , *F^ U and derivatives 
thereof involve an even number of derivatives. Hence, when employing the constant held 
result (H]) for slowly varying inhomogeneous helds, the deviations from the corresponding 
(unknown) exact result are of C^((^-) 2 ). 
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Within the above restrictions, Eq. (JT|) can serve as a starting point to study the effec¬ 
tive interaction between dynamical photons and inhomogeneous backgroimd electromagnetic 
fields. For this purpose it is convenient to decompose the electromagnetic field strength 
tensor F introduced above as F —» F^ u {x) + / MI/ (x) into the field strength tensor of 
the background field F /J,l '(x) = d ,L A u (x) — d u A ll (x) and the photon held strength tensor 


/ M "(x) = d t *a u (x) — d u a tl {x) 12], The effective action can then be compactly written as 


S = S'mw + Sint, where Smw = — 7 / 3^fiv(x) is the Maxwell action of classical electro¬ 
dynamics, with 3 lw (x) = F IJ,u (x) +f ,xu (x) denoting the held strength of both the background 
and the dynamical electromagnetic helds. The additional interaction term S int = 
encodes quantum corrections and vanishes in the limit h —» 0. It can be expanded in terms 
of interactions involving n photons, i.e., ~ f n , with / = f^ v . In particular, to quadratic 


order in / it is given by S int = S^j + + 0(f 3 ), with 


:( 2 ) 


cl 1 ) _ 

D int — 


and = \ ff‘^) gF tg F ^ 


»/ 


I 


X) 


( 2 ) 


The neglected higher-order terms of 0(f 3 ) correspond to interactions involving three or more 
photons. For a somewhat similar approach in a completely different context, cf. 47), |48|. 

Due to the fact that the effective Lagrangian (JTI) is local and thus depends only on a 
single position space coordinate x, all the effective interactions to be induced in the limit of 
slowly varying helds will be local, or - attributing some nonlocality to the spatial derivatives 
- at least “almost local”. 

Equation ([I]) is straightforwardly differentiated with respect to F ^, and the tensor struc¬ 
ture of the resulting expression can be spanned by F^ IU and *F fJjU (cf. also 
convenient to write it as 


491] h We find it 


dC 


d -±F 

dT i,v 


dQ ^ 


( 3 ) 


To keep notations simple, we have omitted any explicit reference to the x dependence of 
the electromagnetic helds. The tensor structure of the second derivative is slightly more 
complicated. It is spanned by six independent tensor structures and reads 


d 2 C 


dF a PdF^ u 


(y9afi9f3v 9cw 9 811 ) 


d£ 

OF 


_* * £ 

7^7 + F^Fafi 77^77 + F/iu Fap-gg^ 

d 2 C 


dC 


d 2 C 

; 8F 2 


{F^*F af 3 + *F )ll/ F OL p ) 


dTdQ 


( 4 ) 


For the reasons given above for the effective Lagrangian, Eq. (|4|) is valid up to corrections 
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Turning to the Fourier representation of the photon field, a p (x) = f k e lkx a p (k) , we write 
the probe photon field strength tensor as f pu (x) = i f k e lkx \k p g ua — k u g flcr ]a a (k). With the 
help of this expression, and taking into account that C is localized, i.e., vanishes at large x, 
such that we can safely interchange the momentum and space integrations, we obtain 

sS=[f(kK(k) and sg = ~ ff a„{k)H<"{k, k')a„(k') , (5) 

Jk z Jk Jk' 

where we have introduced the photon current 


f(k) 


6S int 

Sa^k) 



and the photon polarization tensor in momentum space 


( 6 ) 


n pa (k,k') = 


5 2 S\ 


int 


a =0 


5a p (k)a a (k') 

= {g pp k Q - g pa kP) 


i(k+k')x 


d 2 C 


d F°h j d Fi IU 


{x) 


[k' p g va -k' v g pa ). (7) 


For completeness note that the corresponding polarization tensor in position space, obtained 
by a Fourier transform of Eq. (IT)) , can be represented in an analogous fashion, 


W a (x : x') = 


5 2 S; 


int 


a =0 


8a p {x)a a {x'] 

= (g^d* - <raZ) (g™%, - gre£) 


5{x — x') 


d 2 C 


dF a hdF^ u 


[x] 


( 8 ) 


Tensors of higher rank, describing the effective interaction of n photons can be derived along 
the same lines. 

As the field strength tensor f pv of the probe photon field in momentum space is linear in 
k a , the photon polarization tensor in momentum space is at least quadratic in k a . Taking into 
account that our approach is limited to slowly varying fields from the outset, we have O(k) = 
0{v). Hence, the neglected contributions to the polarization tensor are ~ u 2 C ) ((^) 2 ). 
Employing the shorthand notation ( kF) p = k u F up , ( k*F) p = k u *F Ufl , (, kk') = k p k' p , etc., 
upon insertion of Eqs. fl3j) and (J3J into Eqs. OS]) and (JTJ) we obtain 


f{k) = i 


Akx 


(kFY^Hk’FT | 


( 9 ) 


and 
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n pa (k,k') = 


/ p,i(k+k')x 


dr dr 

(( kk')g pa - k' p k a )— + k'^e^ 


(. kF) p (k'F) a 


d 2 r 

W 2 


+ ( k*F) p {k'*F) a 


dG 

g 2 c 


dG 2 


[( k*F) p (k'F) a + ( kF) p (k'*FY] 


d 2 C 

GTdG 


(10) 


Note that k a j a (k) = 0 and k p U p<T (k, k') = H pa (k, k')k' a = 0, i.e., the Ward identity is 
fulfilled, and correspondingly gauge invariance with respect to gauge transformations of the 
probe photon held is ensured. In the limit of constant background fields, the integrands 
in Eqs. (JTJ) and (TTOj) do not depend on x apart from the overall phase factor e l ( k + k ') x ^ such 
that the integration over x can be performed right away resulting in an overall momentum 
conserving Dirac delta function ~ 8{k + k’). This reproduces the known expression of the 
photon polarization tensor in the constant held and small momentum limit 50]. 

The hve different derivatives of the effective Lagrangian ([I]) with respect to T and G 
appearing in Eqs. (19])- (TTOll can straightforwardly be taken, and their explicit propertime 
integral representations read 


GLl 
GF 

or 

dG 

d 2 r 


a 

2n 

a 

2n 

a 


d 2 L 

dG 2 


d 2 L 

GTdG 


dl 2 2 tt 


ds _ • rr? f ab ascot(bs) , 2l 

s \a 2 + b 2 smh 2 (as) 3 J 

ds _jn£ s f coth(as) cot(6s) 1 6scot(6s) 1 

| 2 ab a 2 + 6 2 sinh 2 (as) + a ^ 1 J’ 

ds 

— e 1 - s 


ab 


abs 2 


(a 2 + b 2 ) 2 \ sinh 2 (as) sin 2 (6s) 


a 

2n 


a 


ds 


-i^-s 


2 as coth(as) 
ab 


i 2 — b 2 


as cot(6s) 
+ 6 2 J sinh 2 (as) 


T (a yy i6) /, 


abs 2 


a 2 + b 2 ) 2 \ sinh 2 (as) sin 2 (6s) 


6 2 


2 a 4 + 6 


A -i 


1 - 2—scoth(as) + _ 

a a 2 a 2 + 6 2 


as cot(6s) 
sinh 2 (as) 


+ (a yy ib) >, 


ds __ iS £ s G 


1 (a 2 - 6 2 )s 2 


(a 2 + 6 2 ) 2 \ 2 sinh 2 (as) sin 2 (6s) 


2 as coth(as) + 


1 (a 2 — 6 2 ) 2 i 6scot(6s) 
6 2 a 2 + 6 2 . sinh 2 (as) 


+ (a yy ib) >. (11) 


Remarkably, for G = 0, i.e., either for orthogonal electric and magnetic helds, or for purely 
electric or magnetic helds, respectively, the propertime integrations in Eq. (TTTj) can even be 
performed explicitly (cf. also [5(, |4l|, l5lj] ). Employing integrations by parts, all integrations 


can be expressed in terms of the elementary integrals (formulae 3.381.4 and 3.551.3 of 52[) 
f^e^' = rr H and / 0 °° f s u e"^ cotli(s) = [2^C(^ f ) - P~ v ]T(v), valid for 


6 




































9?(/3) > 0 and under certain conditions on is, which are rendered irrelevant upon combination 
of these integrals in Eq. (1111) . In result, the nonvanishing derivatives in this limit are 


sc : 

~df 

d 2 C 

W 2 

d 2 C 

W 


s=o 

g=o 

g=o 


^| 4 C , (-1,X) -x[2C'(0,x) -ln(x) + x] - ^ln(x) - ^ j, 
^^{^-x[ 2 C / (°,x) + ln (x) + 2 x( 1 -^(x)) -!] J, 

^ ~~ x ) - ln (x) + x] ~ £ [2-0(x) + X " 1 + !] J, 


( 12 ) 


with 


2e^2|^ h 


for J 7 > 0 
for T < 0 


(13) 


while = Sk\g-o = Here ’ ’I’M = ^lnT(x) is the Digamma function and <(s,x) 

is the Hurwitz zeta function; ('(s,x) = d s ( (s,y). Hence, within the above restrictions, for 
Q = 0 explicit analytical insights into the photon polarization tensor in background fields of 
arbitrary field strengths are possible. 

Another important parameter regime is the limit of generic weak background fields. 
Counting O(a) ~ 0(b) ~ 0(\B\) ~ 0{\E\) ~ O(e) and expanding Eq. (ITT]) in powers 
of < 1 , the propertime integrations can be performed and all prefactors can be de¬ 
termined explicitly. In particular at leading nonvanishing order in such an expansion we 
find 


dC 

dT 

dC 

dg 


a 1 / e \ 2 I 4 T 

7r45vm 2 / j 7 Q 


0 (( 5 ) 4 ). 


(14) 


wherefrom - of course - also J^f, and -^§g can be straightforwardly inferred. Obviously, 
the leading contributions to the photon polarization tensor (ITU]) are of 0((^) 2 ). As our 
approach neglects terms ~ u 2 C>((^;) 2 ), and the vacuum polarization tensor at zero field 
scales as ~ ( k 2 g pa — k p k a ) ~ v 2 0 ((^) 2 ), no zero field contributions show up. 

Finally, we focus on the special case of orthogonal electric and magnetic fields of the same 
amplitude, i.e., E = £ce and B = Ses, with eg • eg = 0, which is of major interest when 
aiming at describing realistic laser backgrounds. In this limit we have T = Q = 0, implying 
that fy 7 = f§ = QJFQQ = 0 - The only nonvanishing coefficients are 


01 = M 2 / 4 

J 7 t45 \m 2 ) 1 7 


(15) 


This significant simplification is due to the fact that for T = Q = 0 the field dependent 
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higher order contributions to the photon polarization tensor scale as ~ m 2 [(-^ I ) 2 0( J ^ I )] n , 
with n G N + (cf. above). The terms with n > 2 are not accounted for when extracting 
the photon polarization tensor from the effective Lagrangian (JT]) , which inherently neglects 
contributions ~ v 2 0((^) 2 J. Hence, in this limit the vacuum current vanishes, j a (k) = 0, 
and the photon polarization tensor is given by 


H pa (k,k') = -y(— 2 

7 T 45 \m z 


j(k+k')x 


4 ( kPyik'Py + 7 ( k*P) p {k'*F) c 


£ 2 


(16) 


where we made use of the definition £F^ = F^ u , i.e., scaled out the local electric field 
amplitude £, rendering the tensor structure in the squared brackets independent of £. In 
turn this tensor structure only depends on the background field alignment to be characterized 
by e ^ = ( 0 ,e B ), = ( 0 ,e B ), = (1 ,eg x e B ), and the probe photons’ in- and outgoing 
wave vectors k' 11 = (a/, k ') and = (ca, k ). 

A particularly convenient representation of the four vectors in Eq. (TT 6 I) is 


{kPy = (kk)e^ - (/ce E )k M , 

(k*FY = (kk)e% - (ke B )k» . (17) 


Assuming the background electromagnetic field to correspond to a plane wave or laser 
field, the unit vector k points in the propagation direction of the laser field. Subsequently, 
we specialize on a globally fixed i.e., the laser’s propagation direction does not vary with 
position and time. Correspondingly, k 11 provides a global reference direction with respect 
to which the probe photon momenta can be unambiguously decomposed into parallel and 
perpendicular components, 

= k^ + k ^_, k^ = (u, k \\), k[[_ = (0, k±) , (18) 

with k\\ = ( k-k)k and kj_ = k — k\\. Let us emphasize that even for fixed the orientations 
of the electric and magnetic fields may still depend on space and time, as it is, e.g., 
the case for circularly polarized electromagnetic fields. 

In order to make the subsequent considerations as transparent as possible, without loss of 
generality we use coordinates where k M = (1, e z ). The directions of the electric and magnetic 
fields can then be parameterized by e B = (cos 0 , sin 0 , 0 ) and e B = e B \ with angle 
parameter 0. For a linearly polarized laser beam we have </> = const., while for circular 
polarization the parameter 0 generically varies as a function of space and time. In these 





FIG. 1: Schematic depiction of a focused Gaussian laser beam (wavelength A) propagating along z. 
The beam focus is at z = 0, where the beam diameter is 2wq. The beam is rotationally symmetric 
about the z axis; r = yhc 2 + y 2 . The Rayleigh range z R is the distance from the focus (along z) for 
which the beam radius is increased by a factor of \/2, i.e., w(±z R ) = \/2wq. The diffraction angle 
is 6 ~ sa. 

Z R 


coordinates, we have 


(kFY = £i(k) cos0 + £ 2 (k) sin0 and (k*FY = K , (19) 

I <p 2 


with £i(k) = (—fc x , k 7 — uj, 0, — k x ) and £%(k) = {—k y , 0, k z — oj, —k y ). Note that these four 
vectors fulfill £]_{k)£ 2 {k) = k£\{k) = k£ 2 (k) = 0 and £\{k)£i(k) = £ 2 {k)£ 2 {k) = ( kk) 2 . 

Using the conventional notations, the electromagnetic field amplitude profile of a linearly 
polarized, focused Gaussian laser pulse (cf. Fig. Q]) in the paraxial approximation (cf. also 
below) can be represented as 


f(r,z,t) 


( z-t ) 2 

£ 0 e h / 2 ) 2 


J^ e -(^y) 2 

w( z) 


cos($(r, z, t)), 


( 20 ) 


with 

<F(r, z, t) = n(z -t) + - arctan(^) + (p 0 , ( 21 ) 

i.e., it generically decomposes into an overall pulse envelope and a modulation term ~ cos(<h). 
The held or polarization orientation, respectively, is fixed by a particular choice of the angle 
parameter 0 = const. 

The adjustable parameters in Eqs. ([ 20 ]) and ()2T]) are the peak held strength £ 0 , the pulse 
duration r, the laser frequency = ^, the beam’s waist size wq in the laser focus at z = 0 
and its Rayleigh range z R = Moreover, w( z) = vj 0 ./ 1 + (^) 2 describes the widening of 
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the beam’s transverse extent as a function of z. The Rayleigh range is the distance from the 
focus (along z) for which the beam radius is increased by a factor of \/2, i.e., w(±zr) = \Z2wq, 
and correspondingly the beam cross section is increased by a factor of two as compared to 
the beam waist. Finally, R(z) — z[l + (^) 2 ] is the radius of curvature of the wavefronts, 
and the term arctan(^) in Eq. (|2TT) accounts for the Gouy phase shift; cf., e.g., 53J. In 
addition we accounted for a phase ipo . The dominant frequency scale governing this type 
of field configuration is Id, which in our approach, counting O(Q) = O(v), is constrained to 
small values of -c 1 (cf. above). 

An alternative representation of Eqs. (1201) and (]2Tj) . which is particularly convenient for 
our purposes, is obtained by elementary manipulations and reads 


£(r,z,t) = £q e 2 ) 2 


w 0 \ 2 1 


w z 


z=±l 


1 _LI IX 




( 22 ) 


The (leading order) paraxial approximation adopted here is valid for small diffraction angles 
9 ~ and neglects terms of 0(9) and It can be systematically improved to 

account for higher order contributions in the parameter 9 54-1561]. In particular, note that 
at leading order in the paraxial approximation, the electric and magnetic field vectors are 
still orthogonal to each other and to the laser propagation direction; cf. also our assumptions 
below Eq. (ITS]) . 

To describe a circularly polarized Gaussian laser pulse we use the following electric and 
magnetic fields, 


(z-t) 


E±(r,z,t) = £ 0 e d/d 2 W ° e [cos($(r, z, t))e x ± sin($(r, z, t))e y ] 


w z 


— So e Em 1 


_(z-t)2 , w N 2 X 

/ u I 


w(z) 


-iZ—)e (1 il ^R ){ ^ )2 e iin ^- t)+ilip0 (e x Ti l ey) (23) 

F_LI 


l = ± 1 


and B± = E± where +/— stand for right-/lefthanded circular polarizations. 

Equation (|22| ) and the second line of (l23l) are particularly suited for an explicit evaluation 
of the photon polarization tensor in high-intensity laser backgrounds. Noteworthily, for these 
field configurations all integrations in Eq. (|T6|) can be performed analytically. Our result for 
the photon polarization tensor in a linearly polarized Gaussian laser pulse reads 


n pa (k,k') = 


45 4 lm 2 ) W ° Zr 2 


4 ( kF) p (k'F) a + 7 ( k*F) p (k'*F) c 


x e -Kl) 2 G+^) 2 


■ p~W 0 (k ± +k ' ± ) 2 

2 g ^ o ( k ±+ k '±) 2 


^w^k ± + k' ± y 
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J ? ^ e -fe«(fc+*')ft e(i ZR (k + £/)£ - >o(*± + *4) 2 ) 

* 2 z=±i 


■ ( 24 ) 


where ©(.) is the Heaviside function, and we used the shorthand notations (k + k')k = 
(k + fc') M «± and (k± + k' L ) 2 = (k± + /d|_) 2 . In this case the electric and magnetic fields point 
in fixed directions and the tensor structures in Eq. (fT6j) are independent of space and time 
and can be factored out. 

The analogous expression for circularly polarized Gaussian laser pulses is 

2 Z^[(fc + fc')ftp g—|t«o(fe_L+fcy ) 2 


X i 22 e [el(k)eW) + 4(k)4W)] 

v W 0\' t -L + kj_) 

- 3 \/f S e-^ )2n[n+Ku+U},)]+il2 *° e ~ l ^ k + k ') k 0 (l ZR (k + k')k - \wl(k± + fcl) 2 ) 

* ^ z=±i 

'[e?(fc)et(*0 -4(k)4(k')\ ± H)'[4W4OT + (25) 


where the subscript ± refers to the results for right-/lefthanded circular polarizations (cf. 
above). As for circular polarization the electric and magnetic held vectors rotate as a function 
of position and time, i.e., 0 —> ±$(r, z, t), in this case we find it more convenient to represent 
the tensor structures in terms of the four vectors e^ and £%. 

Equations (1241) and (J25]) are made up of three generic contributions: One term that is 
independent of the laser frequency 0, and two O dependent terms. The structure of the 
O independent term differs significantly from the others. In the limit of an infinitely long 
pulse duration r, the first term corresponds to an elastic photon scattering process with 
c o = —a/, and the latter terms can be identified with inelastic processes to be associated 
with the absorption (emission) of two laser photons of frequency 0, i.e., co = — (u/ ± 20); 
note that lim r _ ) , 00 |eGi) 2 ^ 2 = y/nS^). 

This generic structure is well known from the photon polarization tensor in plane wave 
backgrounds 42, [43(. In particular at order (^) 2 and for soft electromagnetic fields, the 


plane wave photon polarization tensor comprises the processes of zero and ±20 photon 
absorption/emission from the background field 29]. 


III. CONCLUSIONS 

In this paper we have derived expressions for the photon polarization tensor in inhomo¬ 
geneous electromagnetic background fields, valid for electromagnetic fields which vary on 
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scales much larger than the Compton wavelength and Compton time of the electron, respec¬ 
tively. While for homogeneous backgrounds the photon polarization tensor in momentum 
space can be fully expressed in terms of the momentum transfer (k + k') 1 *, for inhomogeneous 
backgrounds, it generically mediates between two distinct in- and outgoing four momenta k ,fl 
and k G In contrast to previous calculations, we have accounted for the full tensor structure 
spanned by k /fl and AG, rendering our result fully consistent with the Ward identity. 

For most generic electromagnetic field configurations and arbitrary field strengths, the 
resulting expression can be written in terms of a Fourier integral over a single propertime 
integral. In the special case of either orthogonal electric and magnetic fields, or purely 
electric or magnetic fields, the propertime integration can even be performed explicitly. 
Besides, the propertime integrations can also be performed explicitly in the limit of generic 
weak electromagnetic fields. 

A central outcome of our analysis is compact explicit expressions for the polarization 
tensor in realistic laser backgrounds, namely for focused linearly and circularly polarized 
Gaussian laser pulses which are well described by the paraxial approximation, i.e., for small 
diffraction angles, 9 <C 1, and large values for the product of the laser pulse duration and 
frequency rfl 1. Other restrictions for our result to be trustworthy are r c O <C 1 and 
■j— (t c Q ) 2 <C 1, with laser peak intensity I = £ 2 ; recall that the electron Compton time 
(length) is r c = 1.29 • 10” 21 s (A c = 3.86 • 10” 13 m ), and the critical intensity is J cr = (y^) 2 = 
4.68 • 10 29 ^2. All these conditions are typically met for contemporary optical high-intensity 
laser systems delivering multicycle laser pulses. Moreover, the probe photon frequencies have 
to fulfill {t c u>, t c u/} -C 1, and their momenta {A c |/c|, A c |fc'|} <C 1, which is certainly true for 
the whole range from optical to soft X-ray frequency photons available in the laboratory. 

In particular, the latter results are of utmost importance for the investigation of quantum 
vacuum nonlinearities in realistic high-intensity laser experiments, facilitating, e.g., a deep 
and thorough analysis of vacuum birefringence j3| or quantum reflection 281. 
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